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DUST-COOLING-INDUCED FRAGMENTATION OF LOW-METALLICITY CLOUDS 
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ABSTRACT 

Dynamical collapse and fragmentation of low-metallicity cloud cores is studied using three-dimensional 
hydrodynamical calculations, with particular attention devoted whether the cores fragment in the dust-cooling 
phase or not. The cores become elongated in this phase, being unstable to non-spherical perturbation due to 
the sudden temperature decrease. In the metallicity range of 10"^- 10~^Zq, cores with an initial axis ratio > 2 
reach a critical value of the axis ratio (> 30) and fragment into multiple small clumps. This provides a possible 
mechanism to produce low-mass stars in ultra-metal-poor environments. 

Subject headings: cosmology: theory — hydrodynamics — instabilities — stars: formation 



L INTRODUCTION 

It has been suggested that the first stars in the universe were 
very massive, typically > IOOMq (e.g., Abel, Bryan, & Nor- 
man 2000, 2002; Bromm et al. 1999, 2002; Omukai & Palla 
2001, 2003). On the other hand, in present-day local star- 
forming regions, typical stars are low-mass objects < \Mq. 
One of the major factors distinguishing these two modes of 
star formation is the degree of metal enrichment. At some 
point in the history of the universe, metallicity exceeds a crit- 
ical value and the transition to low-mass star formation mode 
is thought to take place. Since the exact value of critical 
metallicity is crucial in modeling such important problems as 
galaxy formation and reionization of the intergalactic medium 
etc., some authors have tried to establish this value. Bromm, 
Ferrara, Coppi & Larson (2001) studied the fragmentation of 
low-metallicity clouds, and concluded that clouds with metal- 
licity > 5 X 1O~''Z0 fragment into smaller pieces owing to the 
fine-structure line cooling of carbon and oxygen. Taking into 
account possible variations in the elemental abundance ratio, 
Bromm & Loeb (2003) derived individual critical abundances 
of C*" and O. Santoro & ShuU (2005) included silicon and iron 
in a similar analysis. However, due to limitations imposed by 
inclusion of only the atomic line cooling, these authors only 
studied the thermal evolution in the low-density regime. The 
mass scale of such fragments is still high > 10- IOOM0, and 
another phase of fragmentation is necessary to produce low- 
mass (< \Mq) objects. Regarding this problem, based on the 
thermal evolution described by Omukai (2000), Schneider et 
al. (2002, 2003) pointed out that another episode of fragmen- 
tation can be caused by dust cooling at higher densities (say, 
> 10'"cm~^). They concluded that dust-induced fragmenta- 
tion can take place even with metallicity of lO'^^'Z© and this 
produces low-mass (< \Mq) fragments. However, their treat- 
ment of dynamical evolution was too simple and their predic- 
tion of fragmentation remained speculative. 

In the previous paper (Omukai et al. 2005), using an up- 
dated thermal evolution and the linear theory for growth of 
elongation (Hanawa & Matsumoto 2000), we evaluated frag- 
mentation mass scale under the assumption that fragmenta- 
tion takes place when the elongation exceeds a critical value 
(~ 1) in the linear theory. However, an approximately spheri- 
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cal core, for which the linear theory is applicable, is not sub- 
ject to fragmentation. Only a highly elongated, filamentary 
(or disk-like) object is gravitationally unstable and easily frag- 
ments (e.g., Larson 1995; Tsuribe & Inutsuka 1999b). Con- 
sequently, a discussion based on linear theory needs to be jus- 
tified by numerical calculations. 

In consideration of the above, here we follow the evolu- 
tion of low-metallicity cores during the dust-cooling phase us- 
ing three-dimensional hydrodynamical calculations and dis- 
cuss the conditions for fragmentation into low-mass clumps. 
Our results confirm that the previous analysis by Omukai et 
al. (2005) is broadly correct: fragmentation does occur in the 
dust-cooling phase and sub-solar mass fragments are indeed 
produced. 

2. NUMERICAL MODELING 

H2-line cooling induces fragmentation of low-metallicity 
(< 1O~''Z0) interstellar clouds around 10'*cm"'', and cloud 
cores are produced. Subsequently, the cores start a so-called 
run-away collapse, where they contract at almost the free- 
fall rate. Here we study evolution of such cores using three- 
dimensional hydrodynamical calculation. Since the cores are 
already self-gravitating, the dark-matter gravity is not consid- 
ered. For simplicity, the cores are assumed to be non-rotating 
and isolated. 

During the run-away collapse phase, a metal-free pre-stellar 
core obeys the self-similar flow with 7 ~ 1.1, as indicated 
by one-dimensional radiation hydrodynamical calculations 
(Omukai & Nishi 1998). Here we choose such a self-similar 
solution as an unperturbed initial state. The initial central 
number density is lO'^cm"-' where the dust cooling is about 
to become effective. 

We add both non-spherical density/velocity and random- 
velocity perturbations to the initial state. As a non-spherical 
density perturbation, the sphere is elongated to a prolate with 
expanding in the z-direction and shrinking in x- and y- direc- 
tions homologously while the central density is kept constant. 
We adopt prolate elongation because we find that the initial 
oblateness leads the amplitude of elongation to saturation in 
the non-linear regime before fragmentation takes place. Elon- 
gation of the core is defined hy £q = a/b-\, where a (b) is 
the short (long, respectively) core axis length. Here, cases 
with an initial elongation of £0 = 0.1,0.32,0.56, and 1.0 are 
considered. Non-spherical velocity perturbation is added to 
the unperturbed flow. The amplitude of non-spherical veloc- 
ity perturbation is assumed to be proportional to the displace- 



2 



Fragmentation of Low-metallicity Clouds 




10 



15 



log n„ 



Fig. 1. — Effective equations of state (EOS) for clouds witli Z/Zq = 
10"'. 10"'' and 10"* (solid). Also shown (dashed) are approximate relations 
(eq. [3 used in the hydrodynamical calculation. 



ment from the sphere and chosen so that the correct linear 
growth (or damping) rate is achieved. Additionally, a random 
velocity field with an amplitude of O.lcj is added, where c, is 
the isothermal sound velocity at the center 

Solving full thermal equations along with hydrodynamics 
is numerically too expensive. Here we take the following 
simplified approach: First, using a one-zone model, we solve 
thermal and chemical processes under the assumption of free- 
fall collapse. The model is the same as in Omukai et al. 
(2005), except the dust model, for which we use dust pro- 
duced by the pair-instability supernova of a metal-free 195M0 
star (Schneider, Salvaterra, & Ferrara 2004; Schneider et al. 
2006). The effective equations of state calculated with this 
model are shown in Figure 1 for metallicities 10"^, 10"^, and 
IQ-'^Zq. Next, for easier analysis and interpretation of the 
result, we approximate the above equations of state using 
the following barotropic relations; the ratio of specific heat 
7 = dlogp/dlogp is given by 



7 = 



1.1 (n<ni) 

0.5 («!<«< na) 

1.1 (n2 < n < rij,) 

1.4 (n3 < n), 



(1) 



where log rti = 6-log Z/Zq, log n2 = ll.l-O.41ogZ/Z0, and 
log n3 = 8.47- 1.21og Z/Zq. In Figure 1, these relations are 
also shown by dotted lines. The first phase (n < «i) corre- 
sponds to the H2 -cooling phase, the second (rti <n < 112) and 
third («2 < « < «3) to the dust-cooling phase before and after 
the gas-dust coupling, and the last («3 < n) to the adiabatic 
contraction after the core becomes optically thick at tij,. 

Godunov-type smoothed particle hydrodynamics (SPH) 
with a second-order Riemann solver are used for hydrody- 
namics and a Tree method is used for gravitational force 
(Tsuribe & Inutsuka 1999a). The number of SPH particles 
is = 1.6 (in Runs A, B, D, F, and G) or 2.3 x 10^ (in Runs 
C and E), which is sufficiently large to meet the resolution 
condition for physical fragmentation (Klein et al 2004). 

3. RESULTS 

Major parameters in this model are the metallicity Z, which 
determines the equation of state, and the initial amplitude of 
elongation £q. The numerical results can be classified into 
two types according to whether the core fragments or not. In 
the following, we show the typical models for these cases. 

3.1. Cores with Fragmentation 



Figure 2 shows the evolution of the core with Z = 10"^ Z© 
and £q= 1 .0. Initially, while the core is close to spherical, 
the evolution of elongation obeys the prediction by the lin- 
ear analysis of a 7 = 1.1 cloud. For central densities n > 
ni(= lO^cm"-'), the temperature decreases rapidly owing to 
the dust cooling, and the core collapses at almost the free- 
fall rate. In this phase, the elongation grows in proportion to 
p™ with m = 0.354-0.5, which is between the rate for spher- 
ical pressure-free collapse (m = 0.354) and that for filament 
collapse (m = 0.5). In other words, the rate m is enhanced 
from the linear value for a sphere owing to a non-linear ef- 
fect. The equation of state begins to stiffen to 7 = 1.1 at 
n = n2i= lO'^'cm"^), where the dust and gas couple thermally, 
and become adiabatic (7 = 1.4) at n = ni{= 3 x lO'^^cm"-'), 
where the core becomes optically thick. However, the core 
continues to elongate owing to inertia until n ^ 10"'cm"-', 
where the maximum value of elongation fmax = 61 is attained. 
A long spindle forms in the high-density core. Around this 
epoch, the spindle stops collapsing in a cylindrical radial di- 
rection and then fragments. Using the scale height of a fila- 
ment// = (2c^/7rG/Jc)'''^, where c, and pc is isothermal sound 
speed and density averaged over the z-axis, the interval be- 
tween each fragment is well described by 2ttH, corresponding 
to maximum growing mode in the linear stability analysis of 
equilibrium filaments (Nagasawa 1987; Inutsuka & Miyama 
1997). The filament fragments into four pieces. Since the 
fragments have an infall velocity along the z-axis, they merge 
with each other at the center, owing to small initial random 
velocities. This is partly because we do not include effects of 
rotation and angular momentum. In spite of the merging of 
fragments at the center, the number of fragments does not de- 
crease, at least within the present calculation, since new frag- 
ments continue to be born on the outskirts. We also studied 
a case with the same parameters (Z, £q), but with a different 
seed for the random number. In this case, positions of density 
peaks were altered, but the distance between each fragment 
remained almost the same. The dynamical evolution was also 
similar. 

In the case of lower metallicity (Z = lO'^Z©, £0= 1 .0), the 
maximum elongation becomes smaller (f max ^ 40) because of 
a weaker effect of dust cooling. Fragmentation is marginally 
observed. There are two fragments. 

3.2. Cores without Fragmentation 

Models with smaller initial elongation exhibit different evo- 
lution. As an example, we show one model with Z = 10"^ Z© 
and £() = 0.32 in Figure 3. In the initial stage (7 = 1.1), 
the evolution of elongation is very similar to examples with 
fragmentation (§3.1). Subsequently, in the cooling phase 
(n ^ "1 = 10"cm"^), the elongation grows. The growth rate, 
however, remains close to the linear growth rate (m = 0.354) 
because of the small initial elongation. Then, the growth 
of elongation slows down at n3 ~ 3 x 10''*cm"^, where the 
equation of state becomes adiabatic. The maximum elonga- 
tion reaches 14 in this case, but this is not enough to cause 
fragmentation. Thereafter, the core becomes more spherical 
with contraction, instead of becoming cylindrical. We also 
followed the evolution in the subsequent accretion stage by 
concentrating on the central region. No fragmentation was 
observed in our calculation. 

In this calculation, the maximum value of axis ratio reaches 
14. This exceeds both the critical value for fragmentation vr 
and maximum growing mode 2tt for fragmentation in the lin- 
ear analysis for uniform infinite filaments. Despite this, the 
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Fig. 2. — Density distribution in the x-y plane on the midplane (top panels), x-z, plane of )■ = (bottom panels) for the case with £o = 1 and Z = 10'^ Zq at 
six different stages. The color scale denotes the density in the logarithmic scale. The maximum number density and time are shown in each panel. 
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Fig. 3. — Same as Fig. 1 but for the case with £q = 0.32 and Z = 10"' Zq. 
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TABLE 1 
Model Parameters and Results 



Run 


log Z/Z0 


So 




log "max 




Fragmentation 


A 


-6 


0.56 


16 


16.8 


1.5 


No 


B 


-6 


1.0 


30 


16.6 


2.7 


Yes 


C 


-5 


0.32 


14 


16.4 


1.6 


No 


D 


-5 


0.56 


27 


16.1 


2.8 


Yes 


E 


-5 


1.0 


61 


16.0 


4.9 


Yes 


F 


-4 


0.32 


34 


15.6 


2.9 


Yes 


G 


-4 


0.56 


56 


15.3 


5.1 


Yes 



NOTE: Model parameters are metallicity logZ/Z© and the initial elongation £q. £max is the maximum value of elongation, which is attained at the number 
density jimax ■ ^^max.L is the maximum value of elongation in the linear theory. 



filament does not fragment because of the central density con- 
centration in our case. Accretion onto the central density peak 
is faster than fragmentation on the outskirts. In the cases with 
£() = 0.1, no fragmentation is observed for any of the metallic- 
ities. 

4. CONCLUSION AND DISCUSSION 

We have studied the fragmentation process of low- 
metallicity (10~^- IQ~'^Zq) cloud cores at high densities (> 
lO'^cm"^). Cores with initial axis ratios greater than about 
two become elongated enough to fragment via dust-cooling, 
even with metallicities of lO"*"- 1O~^Z0. This provides a pos- 
sible mechanism to produce low-mass stars in ultra-metal- 
deficient environments. The results of our calculations are 
summarized in Table 1. As a rule of thumb, fragmenta- 
tion takes place if the core elongation reaches a critical value 
Sail ?J 30. According to linear analysis for a filament of in- 
finite length, namely, without density gradient along the long 
axis, an equilibrium filament is gravitationally unstable and 
subject to fragmentation if the axis ratio is greater than tt (Na- 
gasawa 1987). In our case, the filament is of finite length. The 
density gradient toward the center along the long axis makes 
the critical value for fragmentation several times larger than 
that for a filament of infinite length as discussed in §3.2. In 
Table 1, the maximum values of elongation in the linear anal- 
ysis i^max.L are also presented. We see that if i^max.L exceeds 
^ 3, the elongation of the core will eventually reach > 30 
and the core fragments. It should be noted that the actual 
maximum value of elongation and that in linear theory are not 
reached at the same time because of non-linear effects. The 
linear maximum value of elongation Emax.L is attained at the 



end of the rapid dust-cooling phase (i.e., n = n2), since the 
elongation grows (or decays) if 7 < 1.09 (7 > 1.09, respec- 
tively), in the linear theory (Hanawa & Matsumoto 2000; Lai 
2000). On the other hand, the actual maximum is reached 
slightly after the core becomes adiabatic (i.e., n > n-i). 

Although we have demonstrated by detailed numerical 
work that dust-induced fragmentation indeed occurs under 
some circumstances, some points still remain to be improved. 
To avoid too large a dynamic range, we were forced to take 
the initial state of our calculation at a rather high density 
of 10'"cm"^ and parameterize the elongation of cores at this 
point. The question of how frequently such an initial condi- 
tion is realized remains uncertain. The evolution after frag- 
mentation is also to be studied. The final mass of stars is 
determined as a result of complex dynamical evolution such 
as merging and accretion. In this phase, angular momentum 
would play an important role, although we have neglected ef- 
fects of rotation for simplicity. For example, a disk may frag- 
ment into multiples during the accretion phase (e.g., Saigo, 
Matsumoto, & Umemura 2004). These issues should be in- 
vestigated in a future work. 
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